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The magnetization and hole distribution of Mn clusters in (Ga,Mn)As are investigated by all- 
electron total energy calculations using the projector augmented wave method within the density- 
functional formalism. It is found that the energetically most favorable clusters consist of Mn atoms 
surrounding one center As atom. As the Mn cluster grows the hole band at the Fermi level splits 
increasingly and the hole distribution gets increasingly localized at the center As atom. The hole 
distribution at large distances from the cluster does not depend significantly on the cluster size. 
As a consequence, the spin-flip energy differences of distant clusters are essentially independent of 
the cluster size. The Curie temperature Tc is estimated directly from these spin-flip energies in 
the mean field approximation. When clusters are present estimated Tc values are around 250 K 
independent of Mn concentration whereas for a uniform Mn distribution we estimate a Tc of about 
600 K. 

PACS numbers: 75.50.Pp; 85.75.-d 
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I. INTRODUCTION 



In the diluted magnetic semiconductor (Ga,Mn)As the 
Mn atoms substituting Ga ones act as acceptors that 
simultaneously provide both local magnetic moments 
and spin-polarized p-d-type delocalized holes (accord- 
ing to a simple model one hole per Mn atom) which 
are generally believed to mediate the ferromagnetic cou- 
pling [ISIlSSHHEUHmiil This unic i ue 

feature where magnetism is clearly intertwined with semi- 
conductor properties makes (Ga,Mn)As a particularly at- 
tractive material both for basic research and spintron- 
ics applications. The Curie temperature (Tc) has been 
found to depend dir ectly on the delocalized hole con- 
centration 0, S |H El 111 13 13 El- It has also 
been found that annealing of the molecular beam epitax- 
ially grown (Ga,Mn)As removes Mn interstitial donors 
which increases delocalized hole concentration and con- 
sequently T c [3 13 HI The disorder of the Mn 
atoms either increases Tc further for lower Mn concentra- 
tions [3 13 13 13] or decreases T c [3|. Another mech- 
anism which affects hole concentration (decreasingly) is 
the donation of electrons by the As antisite defects in the 
as-grown (Ga,Mn)As [3 IM 13 13 13 HJ . 

However, an additional mechanism which surely af- 
fects magnetic coupling and the related hole distri- 
bution as well as Tc is the clustering of Mn atoms. 
MnQa-MnQa dimer clusters (denoted henceforth by Mn2 
and called "dimer") and MnQ a - interstitial Mn com- 
plexes have been found in as-grown (Ga,Mn)As samples 
from cross-sectional scanning tunneling microscopy im- 
ages |3l| . Some of these MnQ a - interstitial Mn com- 
plexes may disappear during annealing 20], but com- 
plexes where the interstitial Mn is bound to substitu- 
tional Mn clusters can be very stable |3]. The presence 
of Mn2-clusters is natural because for a characteristic Mn 



concentration of x — 5 % and assuming a random Mn 
distribution the probability that a Mn atom belongs to 
a cluster is 1 — 0.95 12 « 0.46. Structural evolution dur- 
ing post-growth annealing can lead to further cluster- 
ing |l4l Il9l l33j | but no precipitation has been found |33| . 
Recently several first-principles calculations based on 
the density-functional-theory have been performed to 
study the Mn cluste ring and the related changes in mag- 
netism 13 13 Hi HI 13 13 El- The general out- 
come of these calculations may be summarized as fol- 
lows. Substitutional Mn clustering in (Ga,Mn)As is en- 
ergetically favorable up to Mn3- or Mn4-clusters [3E3- 
Large magnetic moments are formed at the Mn clusters 
in (Ga,Mn)As [3 13] [see also Ref. [3] for the case 
(Ga,Mn)N)]. At the same time the p-d hole distribution 
grows significantly at the Mn clusters reducing the rel- 
ative amount of delocalized holes available for the long 
range coupling of the magnetic moments of the Mn clus- 
ters 36, 38]. Interstitial Mn is also found to form a stable 
complex with Mn2 [32] ■ 



The aim of this Paper is to study theoretically the elec- 
tronic and magnetic properties of substitutional Mn clus- 
ters in (Ga,Mn)As paying special attention to the hole 
mediated magnetic coupling between the magnetic mo- 
ments of the Mn clusters using first-principles methods. 
We limit ourselves solely to substitutional Mn clusters; 
interstitial Mn, As antisite or other possible point de- 
fects are not considered although th ey may influence the 
physics of (Ga,Mn)As [3131313121120,121 The Pa- 
per is organized as follows. The computational methods 
are presented in Sec. II, the results for single substitu- 
tional Mn impurity, Mn cluster formation and magnetic 
coupling including discussion are presented in Sec. Ill, 
and the conclusions are drawn in Sec. IV. 



II. COMPUTATIONAL METHODS 

Spin-polarized total energy supercell calculations 
based on the density-functional-theory are performed 
for (Ga,Mn)As. The projector augmented-wave method 
together with the generalized gradient approximation 
(GGA-PW91) for exchange-correlation as implemented 
in the VASP code is employed [H El E2 The pro- 
jector augmented wave method has the advantage that 
it is basically an all-electron method but nevertheless 
almost as fast as the usual plane-wave pseudopotential 
method. However, a slight limitation of the VASP imple- 
mentation is the fact that the core states are kept fixed. 
The projector augmented-wave potentials and the pseu- 
dopotentials for plane waves provided in the VASP pack- 
age were compared with an all-electron full-potential lin- 
earized augmented-plane-wave calculation by performing 
calculations for MnAs (a crystal having the GaAs struc- 
ture but Mn substituting Ga). The projector augmented- 
wave calculation was found to agree closely with the 
linearized augmented-plane-wave calculation but to dif- 
fer both quantitatively and qualitatively from the plane- 
wave pseudopotential calculation [3^ |. 

In our calculations plane-waves up to the 275 eV cut- 
off value are included; total energy in several test sys- 
tems change less than 1 meV as the cut-off value is in- 
creased from 275 to 300 eV or from 275 to 325 eV. Since 
(Ga,Mn)As is a half-metal a dense k-mesh is needed. 
The Brilloin zone integrations are done using the linear 
tetrahedron method with Bloechl corrections [40|]. The 
Monkhorst-Pack (0.14 A -1 ) 3 mesh including the T-point 
was found sufficiently accurate for the Brillouin zone sam- 
pling and is used in this study. One or two Mn clus- 
ters (comprising of up to 5 Mn atoms substituting Ga 
sites) are included in supercells of 64, 96 or 128 atoms, 
corresponding to 2x2x2, 2x2x3 and 2x2x4 cubic zinc- 
blende unit cells, respectively. Thus, the Mn concentra- 
tion in our calculations varies in the experimentally rele- 
vant range [Mn] = 1.6 . . . 7.8 %. Various Mn cluster dis- 
tributions and different spin alignments (e.g. TT vs - Tl) 
are studied. Since relaxation effects were found negligible 
in the 64 atom supercell calculations, only unrelaxed lat- 
tice positions with the fixed experimental lattice constant 
a = 5.65 A are used for larger supercells. The orbital 
decomposition analysis is performed using the standard 
methods implemented in the VASP code ED]- The self 
interaction error and the spin-orbit interactions are ex- 
pected to be small for (Ga,Mn)As [H El E| and 
thus neglected in our calculations. 

III. RESULTS AND DISCUSSION 

A. MnGa impurity in the dilute limit 

The calculated density of states (DOS) for a uni- 
form |46j substitutional Mn distribution in our dilute 
limit of x = 1.6 % (or one substitutional Mn atom in 



the 128 atom supercell) is shown in Fig. ^ for reference. 
The formation of the DOS may be followed using the 
simple branching diagram of Fig. ^ ( a ) ■ The removal of 
a Ga atom creates a vacancy (Vc a ) acting as a shallow 
triple acceptor. The six-fold degenerate p-type i2-level 
of the Vca is occupied by three electrons [471 l48| . When 
the Mn atom (3g? 5 4s 2 ) is placed to the Vc a , the Vc a V 
states and the Mn d states hybridize into the & 
bonding states, the e\ & states, and the i a j & t a ^ 
antibonding states. The three electrons from the Vc a 
and the seven electrons from the Mn atom occupy the 
three bonding states, the two ej states, the three 
i(4 bonding states, and two of the three t a ] antibond- 
ing states leaving one unoccupied state (Fig.^(a)). The 
bands and the DOS are subsequently formed as shown 
in Fig. H(b). The most important feature of the DOS is 
the single hump in the majority spin (f ) channel at the 
Fermi level (Ep) which makes (Mn,Ga)As a half-metal 
(Figs (b) and (c)). As discussed above, the hump is 
formed mainly from the t a -\ antibonding states (Fig. 
(a)). Each Mn atom contributes one hole state in the 
unoccupied part of the hump as well as the net magnetic 
moment of 4/xb [two e-f and two t a ^ electrons, see Fig. ^ 
(a)] (hb is the Bohr magneton). Our calculated DOS 
in Fig. H agrees closely with several other independent 
calculations El El El El El El 

The spin-polarized hole density of the hump consists 
mainly of the Mn d part localized around the Mn atoms 
and the As p part that is delocalized around the As atoms. 
This derealization can be seen in the calculated As p 
projection of the hole distribution shown in Fig. (the 
average number of p holes at the As atoms belonging to 
each coordination shell measured from the closest Mn site 
is given; for comparison the number of holes at the Mn 
site is about 0.13). The peaks at the distances of 2.4 and 
6.2 A belong to the nearest and third nearest As coordi- 
nation shells. Otherwise the p hole distribution appears 
to be quite even. The p hole distribution does not change 
significantly when the Mn concentration is changed (cf. 
the solid and dotted lines in Fig. corresponding to the 
Mn concentrations of 1.6 and 2.1 %, respectively). 



B. Formation of Mn clusters 

The energetically most important calculated Mn clus- 
ter configurations are shown in Fig. [3] The supercells 
used are chosen such that the Mn concentration remains 
the same x — 6.3 %. However, in the case of five sub- 
stitutional Mn atoms in the supercell a slightly higher 
concentration of 7.8 % is allowed due to the computa- 
tional limitations. For the cases of 2-4 Mn atoms in the 
supercell the configuration where the Mn atoms share the 
same neighbouring As site is always found to be energet- 
ically most favorable (see the upper left-hand corners in 
Figs El (a) - (c)). Also, it is energetically favorable to 
form one cluster from two separate components in all 
cases shown in Figs [21 (a) - (c). In the case of five Mn 
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FIG. 1: The calculated density of states (DOS) for the Mn concentration of x = 1.6 %■ (a) Branching diagram for a substitu- 
tional Mn atom. The states of the substitutional Mn atom are formed via the hybridization of the Vgs ti and Mn 3d states, 
(b) DOS. Ef denotes the Fermi level, (c) The magnification of the DOS around Ef- 
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FIG. 2: The calculated As p hole distribution. The average 
number of p holes per As site at each As coordination cell 
is given. The numbers are obtained by integration from the 
orbital decomposition. Note that the orbital decomposition 
as well as the choice of the integration volume are not unam- 
biguous. Nevertheless, the relative magnitudes are directly 
comparable. 



atoms in the supercell the energetically most favorable 
configuration is obtained by attaching the fifth Mn atom 
to the stable tetramer [Fig. 0(d)]. The energies needed to 
separate single clusters into two components (or the bind- 
ing energies for the components) are shown graphically 
in Fig.|2las black horizontal bars which contain also the 
corresponding energy values in the units of meV. In all 
these cases the ferromagnetic order is the stablest mag- 
netic phase. By using the above binding energies we have 
calculated the heats of reaction for the cluster formation 
and found that the optimal cluster sizes are tetramers 
being slightly more favorable than the trimers [3^. Our 
calculation showing that clustering is energetically favor- 
able is in agreement with the experimental result that 
Mn2 clusters appear already in the as-grown (Ga,Mn)As 
samples |3 1| . Further clustering is expected during post- 
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FIG. 3: Calculated Mn cluster configurations and total ener- 
gies. The dark and gray balls denote the Mn and As atoms, 
respectively. Separated components are placed at the maxi- 
mum distance available. All the black horizontal bars show 
the (ferromagnetic) separation (or binding) energies; the cor- 
responding values are given in the units of meV in the black 
bars. The gray horizontal bars with the arrow diagrams show 
the spin-flip energies; the corresponding values are given in 
the units of meV inside the gray bars. All the black and gray 
horizontal bars are drawn in the same scale. 



growth annealing. 



C. Effects of clustering 

The dilute-limit majority-spin hump at the Fermi en- 
ergy Ep in Fig. ^ (c) is found to split when Mn clusters 
are formed, and new narrower unoccupied bands appear 
in the gap. This is shown for the Mn monomer, dimer, 
trimer, and tetramer systems at the constant Mn con- 
centration of x = 6.3% in Fig. 0] At the same time the 
hole density grows at the As atom which is situated in 
the center of the Mn cluster. This was shown for the Mn 
dimer in Ref. j3|| and a similar splitting was also found 
in Ref. The increasing localization at the center As 
atom is reflected in the increasing separation and size 
of the split-off part of the hump in Fig. The As p- 
projection in the split-off part of the hole DOS is seen to 
increase in relation to the Mn d-projection in Fig. 0] and 
even to exceed the Mn d projection in the cases of the 
Mn trimer and tetramer [Figs 0] (c) and (d)] . 

The integrated numbers of the Mn d holes at the clus- 
ter Mn atoms and the As p holes at the center As atom of 
the clusters are given in Table D The results in Tabled 
clearly show that the number of the As p holes grows 
relatively faster than the number of the Mn d holes as 
the cluster size increases. In the case of the pentamer, 
the fifth Mn atom lies outside the tetramer and, in ad- 
dition to the center As atom inside the tetramer, there 
are two equivalent As atoms between the fifth Mn atom 
and the tetramer [Fig. 03(d)]. In this case the holes in 
the tetramer part remain intact whereas the numbers of 
the Mn d and As p holes are seen to grow at the addi- 
tional fifth Mn atom and the two associated As atoms, 
respectively (see the last line of Table H}. We expect that 
an additional sixth Mn atom would cause further hole 
growth at the new center- As-atom. 

The As p component of the spin-polarized hole density 
as a function of the distance to the closest Mn atom is 
shown in Fig. for several different cluster sizes. At 
short distances r < 7 A to the closest Mn atom of the 
cluster, the hole concentration increases with increasing 
number of Mn atoms. However, for larger distances r > 
7 A the hole concentrations approach the same curve. 
Therefore, the long-distance magnetic coupling between 
the Mn clusters depends only on their mutual distance, 
not on the size of the cluster. 

In addition to the different ferromagnetic cluster con- 
figurations above also various other magnetic alignments 
were considered. As mentioned above, the ferromag- 
netic alignment is always energetically most stable. The 
intra-cluster spin-flip energy is relatively high for the Mn 
dimer: 211 meV [3j| and significantly higher inside the 
larger clusters, typically 300 meV and more. The spin 
flip energy of a distant component describes the stabil- 
ity of the long range ferromagnetic state and may also 
be used estimating Tc as will be discussed shortly. The 
calculated spin-flip energies A for the spin-flips of dis- 



TABLE I: Holes inside the Mn cluster. The numbers of the 
Mn d holes at the cluster Mn atoms and the As p holes at 
the As atom which is the center of the Mn cluster are given. 
The numbers are obtained by integration from the orbital 
decomposition. Although the orbital decomposition as well 
as the choice of the integration volume are not unambiguous 
the given average numbers of p holes are directly comparable. 



System 


Mn d holes As p holes 


x (%) 


single Mn 


0.13 


0.03 


6.3 


dimer 


0.17 


0.11 


6.3 


trimer 


0.19 


0.19 


6.3 


tetramer 


0.21 


0.29 


6.3 


pentamer 


0.21 (x4) 


0.28 (xl) 


7.8 




0.19 (xl) 


0.10 (x2) 





tant Mn atoms or entire clusters are given in Table ITU as 
well as in Fig. [3] as gray horizontal bars containing also 
the corresponding energy values in the units of meV. The 
Mn concentrations x in Table ITU are chosen such that the 
closest Mn-Mn distance is constant (10.6 A) except for 
the first row where the distance is 13.8 A. The calculated 
A values vary relatively little, between 57 and 81 meV. 
This variation in the A values exhibits no evident trend; 
e.g. by looking at configurations with one cluster fixed, 
an increase of the second cluster size leads to either in- 
crease or decrease of the A value. These fluctuations may 
be caused by e.g. directional effects [3113 ■ This shows 
clearly that new Mn atoms in the clusters promote the 
holes mainly to the center As atom of the cluster and 
thus do not increase the delocalized hole concentration 
mediating the ferromagnetic cluster-cluster coupling. 
The following expression based on the mean-field ap- 
roximation may be used to roughly estimate Tc [Hsl IHfl . 
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where A is the energy difference between the spinglass 
and ferromagnetic arrangements in the supercell, and N 
denotes the number of the magnetic particles in the su- 
percell. However, we choose to calculate A as the differ- 
ence between the anti-parallel and parallel arrangements 
which is known to be a good approximation |5fj l57j |. 
We treat the Mn clusters as single magnetic particles, 
and calculate A from the cluster-cluster spin-flip ener- 
gies. The number of Mn clusters in the supercell is alwas 
fixed to N = 2. 

First, we consider the simple antiferromagnetic- 
ferromagnetic case. Using Eq.(QJ we estimate that for 
a uniform Mn monomer distribution Tc is 220 K for 
x = 3.1% (Table IHJ and 620 K for x = 6.3% 36]. The 
increase of Tc from 220 to 620 K is mainly due to the 
decrease of the Mn-Mn distance from 13.8 to 9.8 A. If at 
x = 6.3% the Mn atoms were distributed randomly at the 
Ga sites, about 1-0.937 12 = 54% of the Mn atoms would 
belong to clusters. Therefore it is of some interest to 
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FIG. 4: Densities of states (DOSs) around the band gap (upper figures) and the magnifications of the majority spin DOSs 
around the Fermi energy (Ef) (lower figures). Figures (a), (b), (c), and (d) give the DOSs for the Mn monomer, dimer, trimer, 
and tetramer, respectively for the Mn density of x = 6.3%. 
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compare the uniform Mn-dimer distribution at x = 6.3% 
(all Mn atoms belong to dimers) with the uniform Mn 
monomer distribution at x — 6.3%. The dimcrization 
reduces T c from 620 K to 313 K (Table |TTJ) . A similar 
dramatic decrease of Tc is found at x = 12.5% from a 
uniform Mn distribution to the corresponding Mn4 clus- 
ter distribution in Ref. (3^, but at x = 6.3% a slight 



increase in Tc is found from a uniform Mn distribution 
to the corresponding dimer distribution. In Ref. 25] it 
was concluded quite generally that clustering decreases 
Tc, which is in agreement with our findings. 

We follow the same procedure as above to estimate 
Tc from Eq.(^| for further cluster distributions. The 
estimated Tc values are given in Table ITT1 where the Mn 
concentrations x are chosen such that the closest Mn- 
Mn distance is always 10.6 A (except for the Mni+Mni 
case). As noted above A does not depend significantly 
on the Mn concentration and therefore neither docs Tc'- 
the calculated Tc values vary between 220 and 313 K 
(Table [nj. This may be related to the fact that the 
average As p hole distribution at large distances (> 7 A) 
is almost independent of the number of Mn atoms in each 
cluster. Thus, as long as the distance between clusters is 
kept constant, Tc does not vary significantly though the 
number of Mn atoms in the clusters is changed. 

In the mean field approximation, Tc is proportional to 
A/N (Eq.(LH). During cluste ring, A may be expected 
to decrease exponentially with respect to the growing 
cluster-cluster distance while the number of clusters N 
in the denominator decreases much more slowly. Thus, 
as already noted above for the uniform dimer distribu- 
tion, clustering decreases Tc significantly. 

With the present high quality (Ga,Mn)As samples 



TABLE II: Spin-flip total energy differences between two sub- 
stitutional Mn clusters. A is the energy difference between the 
system with one of the two clusters in anti-parallel spin state 
(A) with respect to the system with all spins ferromagneti- 
cally aligned (FM). The 128 atom supercell is used. 



System 


FM 


A 




A (mc 


'V) T c (K) 


x (%) 


Mni+Mni 


T •• T 


T •• 


1 


57 


220 


3.1 


Mn 2 +Mm 


TT •• T 


TT •■ 
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73 


282 


4.7 


Mn 3 +Mni 
Mn 2 +Mn 2 


TTT •• T 
TT ■■ TT 


TTT ■ 
TT - 


• 1 

U 


63 
81 


244 
313 


6.3 
6.3 


Mn 4 +Mni 
Mn 3 +Mn 2 




TTTT 
TTT •■ 


-1 
II 


57 
66 


220 
255 


7.8 
7.8 



where the amount of the harmful interstitial Mn and 
AsQa defects has been minimized one achieves Tc 's in 
the range of 159-173 K |20j,|58|. The mean field approxi- 
mation used here may be expected to give an upper limit 
in estimating Tc- The obtained estimates ranging be- 
tween 220 and 313 K in the case of clustering ( Table ITT|) 
agree quite reasonably with these best experimental val- 
ues. This is consistent with the fact that real samples 
contain clusters [23] which, as was shown above, can re- 
duce Tc strongly. 

IV. CONCLUSION 

Substitutional Mn clustering in the diluted (Ga,Mn)As 
magnetic semiconductor is studied by means of spin- 
polarized all-electron density-functional calculations. 
The clustering is an important factor in the typical Mn 
concentration range from 1 to 10 % because already in 
the case of purely random Mn distribution the probabil- 
ity that a Mn atom belongs to a cluster varies correspond- 
ingly from 0.1 to 0.7. Furthermore, our calculations show 
that cluster formation is energetically favorable. The en- 



ergetically most stable clusters are found to consist of 
Mn atoms that surround symmetrically the center As 
atom. The spin-polarized hole distribution is found to 
get increasingly localized at this center As atom and the 
majority spin DOS hump at the Fermi level to split when 
the cluster size increases from Mni to Mn4. At the same 
time the hole density outside the cluster depends rela- 
tively little on the number of Mn atoms in the cluster, 
especially at larger distances (> 1.2x lattice constant). 
This implies that the long range ferromagnetic coupling 
between two Mn clusters depends relatively little on the 
number of Mn atoms in the clusters. Our calculated 
spin-flip energies confirm this expectation. Using these 
spin-flip energies and the mean field approximation we 
estimate Curie temperatures around 250 K for the vari- 
ous studied cluster distributions independently of the Mn 
concentration. In contrast, we estimate for a uniform Mn 
monomer distribution at x = 6.3 % a high value of 620 K. 
The best present (Ga,Mn)As samples achieve Curie tem- 
peratures of 159-173 K. Therefore, our estimated Curie 
temperature values are consistent with the fact that the 
(Ga,Mn)As samples are hampered by the clustering ef- 
fect. 
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